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Abstract 

The dynamics of two-level systems in an external periodic field are investigated 
in general. The necessary conditions of localization are obtained through analysing 
the time-evolving matrix. It is found that localization is possible if not only is 
the dynamics of the system periodic, but also its period is the same as that of the 
external potential. A model system in a periodic ^-function potential is studied 
thoroughly. 

PACS numbers: 82.90.+j, 03.65.Ge, 33.80.Be 
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Tunneling and localization play a central role in quantum mechanics They are 
the two aspects of one quantum phenomenon. The two concepts are extremely important 
in solid state physics 0. The great achievements of Anderson et al. are focused on the 
localization of electrons while much progress on resonant tunneling has been made since 
Esaki's group proposed the model of superlattice P, It is not until the discovery of the 
coherent destruction of tunneling by Hanggi's group p, |, that attention has been 
paid to the study of localization of a single particle in the double-well potential under a 
periodic acting field. Since this external control of tunneling will be practically valuable 
in many research fields such as laser physics, chemical reaction, macroscopic quantum 
mechanics etc [^, |TU|, |TT], [T^, we need more knowledge about this type of localization 



or the suppression of tunneling. 

To study the localization, the most popular model is a quartic double-well system 
perturbed by a periodic monochromatic field. Such a model has been investigated ex- 
tensively using the Floquet theory |I4| , p!5[] . Because it is not analytically solvable, one 



often recourses to the numerical calculation that sometimes avoids physical insights. Re- 
cently it has been demonstrated p, |T7| that a two-level system shows localization if 
the parameters of the acting field are adjusted, which represents a common feature of the 
double-well system. In a recent paper it is shown that the quantum dynamics of the 



two-level system under a periodic external potential can be mapped to the classical one of 
a charged particle moving in the harmonic oscillator potential plus an magnetic field in a 
plane. The behavior of tunneling and localization is fully described by the radial trajec- 
tory of the particle. This gives us an interesting physical picture although the difficulties 
of mathematical treatments are not lessened at all. 

In this letter, we shall consider the dynamical behavior of two-level systems in the ex- 
ternal periodic field which is antisymmetric with respect to time in one period. Although 
the dynamics has been generally discussed through time-advancing matrix in there 



are still some features unexplored. By the same procedure we will study the evolution of 
the particle in the two states and give the necessary condition of localization. We will 
also derive the explicit results when the external field is a periodic 5-function potential. 

The Hamiltonian of the two-level system in the external periodic potential is 

H = -(Ao/2)(|l >< 1| - |2 >< 2|) + V{t){\l >< 2| + |2 >< 1|), (1) 

where Aq is the energy splitting between the states |1 > and |2 >, V(t) is the coupling 
between them induced by the external periodic driving force, and 

V{T + t) = V{t), V{t) = -V{t + T/2), 
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with T being the period of the external field. We call V(t) a generalized-parity potential. 
Define the left state and the right state as, respectively 

1/ (11 > +12 >)/V2, 



and 

|r (11 > -12 >)/V2. 



The wave vector |\E'(t) > can be expanded in the basis (|/ >, |r >). Denote 

\^{t)>=ci{t)\l>+Cr{t)\r> . (2) 
Then C = {ci{t) , Cr{t))'^ satisfy the equation of motion: 

C = MC, (3) 

where 

_ / -zV{t) zAo/2 \ 
1^ zAo/2 tVit) ) ■ 

Since TrM = 0, the time- advancing mapping or propagator A over a single period is a 



2D area-preserving one (detA=l) [|18]. Apparently this conclusion does not depend on 
the form of V{t) . Defining 

A{t) : C(t) = A(t)C(0) 
for < t < T and A = A(T), we have 

C(nT + t) = A(t)A"C(0). (4) 

Consider three states C(0), C(T/2) and C(T) at three times 0,T/2 and T respec- 
tively. The initial state C(0) is arbitrary except that its components must satisfy the 
normalization condition, |q(0)P + |q(0)P = 1. We can write 



C(T) = ( " f ) C(r/2), (5) 

C(T/2) =(^P ) C(0). (6) 
According to the normalization condition, the matrix elements should satisfy 

\a\' + \P\' = l (7) 
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and 

b\' + \5\' = l. (8) 
Replacing t by —t in ^ and taking its complex conjugate we get 



-iV{-t) iAo/2 
iAo/2 iV{-t) 



C. (9) 



Note that V{t) is antisymmetric with respect to t = T/2. If we move the origin t = to 
t = T/2, then V{t) = —V{t). Inserting it into § and after some analysis one can easily 
find the following relations 

p = a*, q = -7* 
s = -(3\ t = 5*. 



Thus, we have 



( \a\'-m' (36* -aY \ ...^ 



Together with ^ and ^, we obtain 

A=(,r -b + ^c\ 

where a, b and c are real numbers determined by the system and + 6^ + = 1 from 
detA=l. 

Because —2 < TrA = 2a < 2, the dynamics of the system is strongly stable [P^. 
In other words, the future behavior of the two-level model is insensitive to the initial 
condition. Roughly speaking, localization means that the state of the system evolves 
near the initial state. To investigate the rule of localization, thus, we can choose any 
initial state as we like. But what we are interested in practice is related to the evolution 
of either the left state or the right one. Therefore, we suppose C(0) = (1,0)^ in the 
following discussion. The state at the nth stroboscopic point nT is given by 

C{nT) = A"C(0). (11) 

The power of A can be determined by the Caley-Hamilton theorem (see, e.g. |19||): 

A"-i = P„_2(a)A - P„„3(a)I, (12) 

where P is the Chebyshev polynomial and I the identity matrix. Denote a = arccosa, 
then Pn = sin[(n + l)cr]/ sin cr. 

After a simple calculation we obtain the modulus of Cr{nT): 

\cr{nT) \ = \sm{na) \ . (13) 
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Localization requires that for any integer n the probabihty of finding the system in the 
right state should be small, i.e., \cr(t)\ « 1. This is only a qualitative statement. Quan- 
titatively, we define localization as \cr\ < 1/2 for all the time. As a necessary condition 



for localization, therefore, we must have \cr{nT)\ 



sin(ncr)| < 1/2. Mathematically if 

(14) 



where P and Q are two mutually prime integers, then the set {sin(?T,cr) : n^N} is dense 
in the interval [0, 1]. In other words, for arbitrary a there must exist such an integer n 
that sin(n(T) is very near to 1, as a consequence, localization is destroyed. To guarantee 
localization, thus it is needed 

for some integers P and Q, i.e., there must be a set of integers m to satisfy 



sm ma 



0. 



Suppose m > 1 and m is the smallest one in the set. Thus a can be written as 

cr = — , 

m 

where / and m are mutually prime integers. If m is even, then 



(16) 



(17) 



sm 



which is forbidden for localization because mT/2 is also a stroboscopic point. If m is odd, 
since m and / is mutually prime, there is an integer n less than m such that mod(n/, m) = 



Then we have 



sm mcr 



sm 



(m — 1)tt 



1 

> -. 
2 



(19) 

1 localization may arise. 



Taking account of |I8| and |I9|, we conclude that only if m 
Namely, we have proved that 

Theorem: The necessary conditions for localization of two-level systems in external peri- 
odic fields are 1) its dynamics should be periodic; 2) the period should equal to that of 
the external fields. In mathematical language, these conditions are simply described as 



A = ±1. 



(20) 



This result strenghthens the conclusion in reference [0 in which only the first local- 
ization condition is given. Note that this observation was obtained by GroBmann and 
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Hanggi in |p, using the Floquet approach. In the following, this theorem will be applied 
to the treatment of the periodic 5-function potential or pulse acting on a two-level system. 



The 5-function potential is very important in physics because it is simply representing 
the most physical features of other acting potentials. We find that it is difficult to solve ^ 
directly, and so we treat the 5-function field as the limited case, e — > 0, of the rectangular 
potential or pulse, namely. 



Vit) 



if < t < T/4 - e 

Vo/2e if T/4 - e < t < T/4 + e 

if T/4 + e < t < 3T/4 - e 

-K)/2e if 3T/4 - e < t < 3r/4 + e 

if 3T/4 + e < t < T 



and V{T + t) = V{t) (see Fig. 1). If we define the time-advancing matrix T of one 
(5-function potential Vo6(t) as 

C(0+) =TC(0-), (21) 



we can derive 



r 



exp{iVo) 
exp(— iVo) 



f22l 



The time-advancing mapping A turns out to be 

1 - 2 sin^ ^ cos^ Vo - sin 2Vo sin 



A 



sin2Vb sin 



^An + JsinlAllcos^V^n 



^ + i sin ^ cos2 Vo 



l-2sin2^cos2Vo 



From the definition of a we have 



cos cr 



1 



. 2 ^'^O 2 

2 sm — ■ — cos Vn. 



Using the theorem above, we find the necessary condition for localization as 



sm 



cos Vn = 0, 



or 



sin Vn = 0, 



cos ■ 



tA, 



0. 



(23) 
(24) 

(25) 
(26) 



If this condition is satisfied, and if the system is localized in one single period, then 
localization will dominate the dynamics for the whole time. In order to find the sufficient 
condition for localization, therefore, we only need study the dynamical features over one 
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period. To this end, it is required that the modulus of q be never less than a half in one 
period. This requirement imposes constraints on Vq and T respectively, i.e.. 



/ Vo = {n + i)7r, 

\ < TAo < 2n, ^ ' 

where n is an integer. This gives the sufficient conditions of localization for a two-level 
system in the periodic 5-function field. 



Returning to |2^, we see clearly that the applying 5-function pulse only brings about 
changes in the phases of the left and the right states. The former is increased by Vq and 
the latter decreased by Vq. Since the evolution is strongly dependent on the difference 
of the two phases, we can always choose a series of appropriate 5-function potentials to 
obtain localization. For example, let us apply a periodic potential consisting of a pair of 
(5-function pulses with opposite signs in one period to a two-level system. The two pulses 
are assumed to take effect at one fourth period and three fourths period (see Fig. 2). 
Suppose the system is in the left state at the beginning. As shown above, if and only if 
the positive and the negative pulses have the same amplitudes, |Vo| = {n + |)7r and the 
period T is not larger than 27r/Ao, then the system is localized and the probability for 
the system to stay in the left state will not be less than pi^min = cos^(TAo/8). Noticing 
that T < 27r/Ao, we find that pi^mm is a monotonically decreasing function of T, which 
indicates that the smaller T is, the larger pi^min- In other words, we can realize strong 
localization by taking short-period 5-function pulses. The varying of the probability ac- 
cording to time for different periods is illustrated in Fig. 3. It should be stressed that 
although the period of the dynamics is identical with that of the external field, the period 
of the evolution of the probability in the left state is a half. In fact, this is a common 
feature of the two-level system perturbed by any external potential possessing the gener- 
alized parity, which can be shown via symmetry analysis . 



In conclusion, we have derived the necessary conditions for a two-level system driven 
by a periodic field to be localized: not only must the system evolve periodically, but also 
the period is identical with that of the external potential. The necessary conditions can 
always be satisfied by adjusting the period and the strength of the applying field. Under 
these conditions the system will get localized if and only if it is localized over one single 
period, which becomes possible through modifying the parameters of the external field. 
In the case of 5-function pulses, the necessary and the sufficient conditions for localization 
are acquired analytically. In principle, any applying potential can be imitated approxi- 
mately in terms of a series of 5-function pulses, and so the general localization conditions 
can be investigated on the basis of the results from our 5-function model. Besides, we 
expect that this simple theory can be employed to study the process of enantiomerization, 
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which is of much significance in understanding the origin of hfe. 
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Captions of figures 

Fig. 1 An applied field consisting of rectangular pulses. It becomes the 5-function 

field by taking the limit e ^ 0. 

Fig. 2 Schematic graph of applied 5-function pulses. 

Fig. 3 Calculated probabilities in the left state versus time for three periods. Here 
the solid line corresponds to TAq — Air/ 5, the dot dashed line TAq — 4n/3 and the 
dashed line TAq = 2tt. 
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V{t) V{t) 

V(t) V(t) 

time (unit: 2/Ao) 2T T/2 

T/4 3T/4 2T 2e VQ/2e 

probability in the left state V(t) —Vq 

time (unit: 2/Ao) T T/2 

T/4 3T/4 T 2e V2e 
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